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ABSTRACT 


The use of thermoluminescent dosimeters (CaFg 
and LiF) for the determination of gamma-heating m 
reactor materials has been assessed m the light of tem- 
perature measurements across certain specially designed 

60 

experimental assemblies. A *+0 KCi Co source was used 
for irradiation. The transport-theory code, ANISN, was 
employed for estimating gamma-spectra within the assem- 
blies to enable the application of general cavity- 
lonization theory for the deduction of gamma-heating 
rates from TLD measurements. The measured heat-generation 
rates were fed as input to a three-dimensional heat con- 
duction code, HEATING, to obtain various temperature 
drops across each assembly. These were compared with 
actual thermocouple measurements, reasonable agreement 
being obtained m each case. 
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CHAPTER 1 
INTRODUCTION 

In the design of nuclear jpower reac tor s, 
knowledge of gamma-ray heating effects is needed both 
for determining the cooling requirements for the reactor 
components and for predicting where thermal stresses 
will be important. If core, blanket and control-rod 
cooling requirements are to be correctly assessed ard 
optimum shielding specified m the design of a fast 
reactor, it is necessary to be able to determine the 
amount of energy that will be deposited by gamma-ray 
interactions in the various regions Gamma-ray heating 
measurements m zero-power critical s allow one to deter- 
mine the adequacy of calculational techniques and gamma- 
ray-production cross-section libraries. They also provide 
gamma- ray heating profiles for a particular design and 
for regions not amenable to calculation. 

A number of devices have been developed for 
gamma- ray dosimetry, e.g. cavity-ionization chambers, 
scintillation detectors, chemical dosimeters, solid- 
state integrating dosimeters, photographic emulsions 
and calorimeters. However, the particular environment 
of a reactor and other operational considerations 
limit the choice that one has. Measurements have to 
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be made m a mixed neutron/gamma radiation field in 
moderate fluences and m a variety of materials. It is 
also often necessary to perform simultaneous measurements 
m a number of different locations with probes that intro- 
duce minimum amount of perturbation. The thermolumines- 
cent dosimeter (TLD) has proved to be the most suitable 
choice for gamma-ray heating measurements. TLD materials 
and techniques have developed over the last decade or so 
to an extent whereby a variety of dosimeters have become 
available for application to a wide range of exposure 
conditions . 

There are two special areas of concern in the 
application of TLD techniques to the measurement of gamma- 
ray heating in reactors, viz, (i) the relationship between 
the energy deposited by gammas m the dosimeter and the 
heating m the surrounding medium, and (11) the neutron 
sensitivity of the TLD material. The first effect is a 
function of the composition of the TLD and the surrounding 
material, the size of the TLD and the gamma-ray energy 
spectrum at the location under consideration. Since the 
gamma spectrum is not generally known accurately and 
since TLD materials are usually quite different from the 
materials m which gamma-heating measurements are to 
be made, it is important to choose the TLD suitably so 
as to minimize the uncertainty. As regards the second 
aspect, all TLD materials exhibit some response to 
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neutrons since ion recoils from elastic scattering, as 
well as gammas and/or other reaction-products from 
neutron capture, can cause ionization m the TLD material . 
Neutron-induced contributions to the TLD response must 
therefore be suitably corrected for in experimental gamma- 
heating studies m reactors. 

In the present work, the adequacy of gamma- 
heating measurements made with CaP^ and L;lF TLD's has been 
directly assessed m the light of temperature measurements 
made using thermocouples. Special experimental assemblies 
were designed to enable such an assessment-steel, aluminium 
lead and teak wood being the materials considered. A pure 
gamma field provided by a ^-0 KCi cobalt-60 source was 
employed for the present study, the energy absorbed by the 
TLD's being related to the energy deposited m the surround- 
ing material through the application of general cavity 
ionization theory. The latter required an estimation of 
the gamma- energy spectra at different location within 
each experimental assembly, and these were obtained from 
multigroup one-dimensional photon transport calculations 
using the AWISN code. With gamma-heating rates at 
different locations m a given assembly thus estimated 
from TLD measurements, a three-dimensional heat-conduction 
calculation was carried out using the code HEATING to 
obtain steady-state temperature drops across certain 
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specific locations. These calculated temperature drops 
were compared with actual thermocouple measurements made 
by irradiating the assembly for a period long enough to 
achieve steady-state temperatures. 

Chapter 2 gives a general review of TLD prin- 
ciples and earlier applications. The present approach 
and design of experiments is discussed m Chapter 3 5 
which also describes the analytical methods used, viz. 
the ANISN and HEATING codes. Particulars of the various 
experimental facilities and equipment used m the present 
study are given m Chapter If. Results are presented and 
discussed in Chapter 5> while the final Chapter summari- 
zes the conclusions and indicates the scope for further 


work , 



CHAPTER 2 


GENERAL REVIEW 

Thermoluminescent dosimeters (TLD’s) provide a 
useful means for measuring gamma- ray heating m critical 
facilities and shielding materials because of their high 
sensitivity, small size? easy handling, broad dose-measure- 
ment range, low neutron sensitivity, advanced state of 
development and adequate reproducibility/ accuracy. Thu s 
chapter gives a general review of TLD applications in 
reactor engineering. 

2.1 TLD PRINCIPLES 

The basic principle of the TLD is quite simple. 
The thermoluminescent (TL) material contains a number of 
electron or hole traps, stable at room temperature. An 
ionizing particle in passing through the material raises 
electrons into the conduction band and some of these 
electrons, instead of falling directly back to the ground 
state, are captured in trapping centres at some inter- 
mediate energy level. A trapped electron cannot get to 
a lower energy state except by first returning to the 
conduction band. When the material is heated to appro- 
ximately 300 °0, the electron (or hole) is released 
from the trapping centre by thermal ionization and 
moves through the crystal until electron-hole annihi- 
lation takes place at a recombination centre. During 
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this process thermoluminescent emission occurs which can 
be detected by a photomultiplier tube. The number of 
photons released, and thus the intensity, are related to 
the material’s exposure. Figure 2.1 (a) depicts the 
trapping of an electron/hole during exposure, while 
Fig. 2.1 (b) depicts the emission of the TL photon on 
heating . 

2.2 SUMMARY OF EARLIER APPLICATIONS 

Extensive developments m TL materials and 
techniques during the last ten years have resulted in 
dosimeters that are applicable to a wide range of exposure 
conditions and are being increasingly employed for per- 
sonnel dosimetry (1, 2) . TLD techniques have also been 
applied to reactor studies (as discussed below) and to 
other mixed radiation environments such as those encoun- 
tered m weapons testing. For gamma-heating measurements 
with CaF 2 • Mn dosimeters m various reactor materials, 
the dose/ response relationship may he determined with an 
accuracy of typically 3$, while for LiF errors are gene- 
rally larger because of the greater energy dependence. 

Calorimetric measurements of gamma-ray heating 
rates m critical assemblies are very difficult, since 
the gamma flux at a typical power level of 50 - 100 
watts would cause a temperature rise of only about 
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o 

10 C/hr, assuming perfect insulation of the calorimeter. 
Consequently, the possibility of indirect measurements 
was considered, and it was concluded that an ionization 
detector would provide the desired information with the 
requisite sensitivity. The TLD was chosen, m preference 
to ionization chambers, because of the inherently smaller 
probe size and the possibility of making simultaneous 
measurements at a number of different positions. Photo- 
graphic emulsions were also considered m the early 
stages, but TLD's were soon recognised as being consider- 
ably more advantageous. 

Stanford and Johnson (3) used CaF^ * hn powder 
to measure gamma-ray heating in a thermal critical faci- 
lity. Measurements were made in the vicinity of the 
core at operating powers of 20 - 300 watts. The Bragg - 
Gray principle of cavity ionization (Sec. 2.3) was used 
m analyzing the data, e e rime nt ally determined correc- 
tions being applied to account for the gamma-ray spectral 
effect and neutron sensitivity. The reproducibility of 
the measurements was typically about 5 - 10$, Absolute 
results were subject to an additional uncertainty of 
+ 15 % m the determination of the reactor power. 

Adamson, et. al. (4) used LiF powder to 
measure gamma-ray heating as a function of position 
within a single lattice cell at the centre of ZEBBA-6, 



a zero-power fast reactor assembly, BouLette and. Bunch 
( 5 ) used. LiF ribbons (3 x 3 x 0,9 mm) for checking the 
shield, design of the Fast Flux Test Reactor (FFTR) . In 
further support of the FFTR programme, Stanford, et al 
(6) exposed^LiF rods (1 x 1 x 6 mm) in the ZPR-9 critical 
assembly, LiF TLD's were also used for c!os^ mappings, 
single-cell dose heterogeneity and blanket-reflector and 
reflector-core interface studies for the Experimental 
Breeder Reactor (EBR - II) series of ZPR - 3 critical 
assemblies (7), 

Absolute gamma-ray heating was investigated by 
Simons and Olson (8) m Assembly-2 of the Zero Power 
Plutonium Reactor (ZPPR) . The first series of experi- 
ments was devoted to obtaining axial and radial absorbed- 
dose maps for the LMFBR Demonstration Benchmark Critical. 
The second and third series of experiments determined the 
heating measurements inside LMFBR control-rod assemblies. 
All measurements were made witb/LiF TLD 1 $ (enriched to 
99*993^ ^Li) encapsulated m stainless steel. Fast- 
neutron spectra were measured inside the plutonium - 
fuelled critical assembly near the centre of the inner 
and outer core zones by proton-recoil proportional 
counters, and also calculated by the DOT code. (All 
neutron cross-sections needed for the DOT code were 

p 

obtained from library tapes prepared using MC , CALHET 
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and CCXS) . The gamma-ray source m a given photon group 

due to neutron induced gamma- ray-pro duct ion reactions 

was calculated hy the POPOPt- code. Finally, using 

Compton - effect, photoelectric and pair production cross- 

sections from a MUG library tape, absorbed - gamma 1 doses 

were calculated m rads/sec from the relation 

N 

D = 0 

3=1 

o 

where C is the conversion factor 1 .602 x 10” gm/Mev 

0 is the photon flux in the 3 th group (3 = 1 , * 

J 

2 

m photo ns/ cm - soc 

(E ) is the average photon energy of the 3 th group, 
av 3 

in Mev 


("yfvO is the mass energy absorption coefficient of 

' 2 

the material for gamma-energy E m cm /gm. 

Stv * 

Excellent agreement was obtained between calculated 

n 

absorbed doses and results of measurements with LiF TLD r s< 


It should be mentioned that a few calorimetric 
measurements of gamma - ray heating in fast reactors have 
also been recently reported, e.g», by Price, Laskiewicz 
and Lowery (9) and by Reilly and Peters(IO) . Two 
different types of calorimeters were used - an active 
calorimeter using thermocouples, and a passive calori- 
meter termed as the gamma expansion difference monitor 
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(GEEM). Hie latter system employed thermal expansion 
difference (TED) monitors, each consisting of a small 
cylindrical stainless steel Capsule completely filled 
with sodium. 

Another type of TLD application m reactor 

engineering, i.e. apart from gamma-heating measurements, 

has been the determination of m-core thermal-neutron 

fluxes (11) , LiF crystals enriched with either ^Li or 
n 

Li allowing the application of two alternative neutron/ 

gamma discrimination techniques. In the first method, 

both types of crystals were employed. The ^LiF (95*62^ 

6 

Li) dosimeters, due to the large thermal neutron cross- 

6 

section for Li, detected both neutrons and gamma radia- 
tion, while ^LiF (99*993% ^Li) dosimeters, possessing 
negligible thermal -neutron attenuation characteristics, 
monitored the gamma component only. The dosimeters were 
inserted into the reactor for a known time and read on a 
commercially available detection system, the difference 
m the response of the two types of dosimeters yielding 
a direct measure of neutron flux. The second technique 
used bare and cadmium-covered ^LiF dosimeters. The bare 
crystals detected reactor gammas, while those encapsula- 
ted m cadmium measured reactor gammas plus capture 
gammas from the ^Gd (n, y ) reaction. The difference 
m response provided the capture-gamma contribution, 
which was proportional to the neutron flux. 
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2.3 CAVITY-IONIZATION theory 

When a dosimeter, e.g. an ionization chamber or 
TLD, is introduced into a medium for measurement of dose, 
the energy absorbed per unit mass m the medium is usually 
different from that absorbed by the dosimeter since the 
dosimeter and surrounding medium are generally different 
materials. Cavity ionization theory (12) provides the 
relationship between the energy absorbed per unit mass m 
the medium, E m , and that m the dosimeter, E^, m the form 



where f is, xn general, a function of the (a) gamma ray- 
energy, (b) composition of the surrounding medium, 

(c) composition of the TLD material and (d) size of the 
dosimeter. 

A simple theory of cavity ionization results if 
one assumes the followings 

(l) The primary electrons created by photon inter- 
actions lose their energy continuously. 

(n) The cavity does not perturb the equilibrium 
electron energy distribution established m 
the medium, and 

(ni) Photon interactions in the cavity material 
are negligible. 



The latter two assumptions imply that the size of the 
cavity is small compared to the electron ranges m the 
cavity material. If one primary electron is generated 
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per unit mass of the medium with energy T q , and if 
electron equilibrium exists, the energy absorbed per unit 
mass of the medium is T . The energy deposited m the 
cavity is 
T 

j[ T)<^i) D dI (2.2) 


where (T , T) is the electron slowing-down energy 

distribution characteristic of the medium and, according 

to assumption (i), is the reciprocal of the mass stopping 

1 dT 

power of the medium. ( )p is the mass stopping power 

for the cavity material. The energy dissipated per unit 
mass of the dosimeter material divided by the energy 
dissipated per unit mass of the medium is given by 


f 


s 


(T ) 
o 



( 1 dT 
_ •£ dx 


D 


( 


1 

f 


dT 

dx 


dT 


m 


( 2 . 3 ) 


If the energy dependence of the mass stopping powers is 
neglected, the above equation reduces bo the well known 
Bragg - Gray relationship. To determine f (T^ ), the 
f _ value associated with a giv^n monoenergetic gamma- 

S j f 

ray of energy T^ , - it a^ragmg over the three 



primary electron spectrum components produced by the 
gamma-ray is carried out, giving 
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f s (T r ) 


F (T v ) f + F 
p a p s c 



f + F 
c s pp 




c 


( 2 . 4 ) 


where the subscripts p, c and pp denote photoelectric, 
Compton and pair production interactions respectively, ^ 
is the interaction probability and f , the corresponding 

P 

f term, 
s 


The small-cavity assumptions made above are not 
valid for typically sized solid-state dosimeters exposed 
m reactor gamma-ray spectra. The electron spectrum m 
the cavity is not identical to that in the medium, and the 
effect of gamma-ray interactions with cavity material is 
not negligible. In the limit of a large cavity (the 
range of the most energetic electrons is small compared to 
the cavity dimensions), the energy deposited m the cavity 
is simply 





-&en \ 
f 


( 2 .?) 


where N = 

h = 

/ >&en \ 

1 f ; D“ 


gamma- ray fluence 
gamma- ray energy 

the mass absorption coefficient for the 


dosimeter material 



A similar expression holds for the surrounding medium. 
Thus in the large-cavity limit, the energy absorbed per 
unit mass of the dosimeter material divided by the energy 
absorbed per unit mass of the medium is given by 

k r 'd 


h < T * > 


/ Men v 
^ f ; m 


( 2 . 6 ) 


A comprehensive theory of cavity ionization must 
reduce to the small-cavity result at large electron ranges 
and to the large-cavity result at small electron ranges. 
Burlm ( 13 ) proposed an approximate method of solving the 
general problem by considering that the presence of the 
cavity in the medium modifies the electron spectrum m two 
ways, viz. 


(1) the electron spectrum generated by the photons 
m the medium and emitted from bhe wall of the 
cavity will be partially or completely absorbed 
as it crosses the cavity and 
(11) the photons will generate 111 the cavity an 
electron spectrum which will build up and 
possibly reach equilibrium at increasing dis- 
tances from the cavity wall. 


To a first approximation, it may be said that the elec- 
trons are attenuated exponentially without any change in 
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the spectral distribution, the effective mass absorption 
coefficients being assumed the same for the "wall spectrum" 
and the "cavity spectrum". Thus the electron spectrum in 
the cavity at a distance x from the wall is reduced by 
the factor e _/ ^ x , and the electron spectrum generated m 
the cavity will have risen to (1 - of its equili- 

brium value • ft is the effective mass absorption coefficient 
and for the case of interacting gamma- rays of energy T^ , 
is given by 


/3 = — 3— nyr cm 2 / gm (2.7) 

(Ty ~ 0.036) 1 ’ 4 - 0 

g, the average path length of electrons crossing the 
cavity, is given by 


g = 


4 V 


( 2 . 8 ) 


where V is the volume and S is the surface area of the 
dosimeter. Then, on average, bhe wall spectrum will be 
reduced by a factor 


f 


,-/6 x 


dx 


r 


dx 


i - dEl 
/?6 


= a 


(2.9) 


The cavity spectrum will build up to a fraction of its 
equilibrium value given by 

r g 

J (1-e“^ x )dx r -/3g 

2 - t e/ . ° - 1 = (1 _ a) (2.10) 


f 


dx 


/3 g 
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The energy absorbed, per unit mass of the desimeter 
material divided, by the energy absorbed per unit mass of 
the surrounding medium may finally be approximated by 

f (Ty ) = d f s (Ty ) + (1 - d) f ± (Ty ) (2.11) 

Equation (2.11), the general-cavity theory result, reduces 

to equation (2.3) for the case of a small cavity 

(g~~ >-0, d = 1) and to equation (2.6) for a large cavity 

(g-^ co ) d = 0). This general result forms the basis for 

carrying out gamma-heating measurements. For an incident 

gamma spectrum described over N energy groups, i.e. 

-1 

J3 , i = 1,2, ... N, the effective jr- factor to be used 
for converting the TLD dose to energy absorption m the 
medium would be 

Sum (~ * 0 ) 

1 - ( 2 . 12 ) 

* Sum( 0.) 

i 

where f is evaluated at the mean gamma energy for 
the ith group. 

2.4 FACTORS EFFECTING RESPONSE-DOSE RELATIONSHIP FOR TLD * S 

A TLD is not an absolute radiation detector. 

It is thus necessary to determine empirically the rela- 
tionship between the TL response and the energy 
deposited m the dosimeter by exposures m known fields 
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of radiation. The relationship between bhe TL response 
and E d may be a function of several variables such as bhe 
TLD material, the dose, the dose rate, the energy of the 
gamma-rays which produced the dose, the heating rate, 
the time between the irradiation and readout, environmen- 
tal factors before, during and after the irradiation, 
and the pre- and post- annealing history of the dosimeter. 

The dose-response curve of CaF^ : Mn is linear 

up to the point at which saturation occurs . LiF exhibits 

■3 

a sup r alme ar behaviour beginning at about ICu rads, below 
which the response is linear (2). 

The energy response is the dependence of the 
dose/response curve on the energy of the gamma- rays 
which produce a given dose. For GaFg * Mn, the response 
per rad appears to be independent of the gamma-ray energy, 
whereas for LiF there are certain finite radiation - 
quality dependent effects m the TL response at low 
electron energies. 

Both LiF and CaFg • Mn show evidence of fading. 
The amount of fading is less than 10% per month and 
appears to follow an exponential behaviour. If all the 
dosimeters used m a particular experiment are always 
read out at a fixed time after irradiation, corrections 
for fading effects become unnecessary. 
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The LiF response is somewhat more sensitive to 
the dosimeter's pasb history than is the CaF^ : Mn response. 
The response of LiF significantly decreases if the phosphor 

If 

has been exposed previously to radiation greater than 10 
rads. This difficulty can be overcome if the calibration 
of the dosimeter is carried out from the same batch* With 
commercially available readout devices and dosimeters, it 
is possible to obtain 2 to 3 % reproducibility for a given 
dosimeter, as well as, from dosimeter to dosimeter within 
the same batch. 

2.5 NEUTRON SENSITIVITY 

When carrying out gamma-heating measurements m 
a reactor environment, the sensitivity of the TLD response 
to fast and thermal neutrons has to be assessed and correc- 
ted for if necessary. 

Thermal neutrons create ionization m the mate- 
rial either through charged-particle reactions, such as 

^Li (n, c< ) % m a LiF dosimeter, or through gamma- 

26 

producing reactions, such as ^Mn (n, y ) Mn m 
CaFg * Mn dosimeters. For gamma- ray heating determination 
m a thermal reactor environment, the CaF 2 • Mn dosimeter 
is preferable because it shows minimal sensitivity to 
neutrons ( 14 - ). The sensitivity of various TLD’s bo 
thermal neutrons has been reported by Busuoli and 
Cavallim 05) and is reproduced below; 
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Material 

ID 

Response/10 cm 

response/R of 

LiF (TLD-600) 

1930 

LiF (TLD- 100) 

350 

Lig Oy % Mn 

310 

LiF (TLD-700) 

1.5 

CaF 2 ; Mn 

0.6 

CaSO^ : Dy 

0.5 

BoO 

0.3 


(relative to 
Co gamma- rays) 


TLD's exhibit a thermoluminescent response to 
fast neutrons. The elastic scattering of fast neutrons 
by the lattice ions produces ion recoils. The slowing 
down of these recoils creates electron-hole pairs, which 
may then become trapped. The slowing down also leads to 
the production of defects - ions are displaced, and m a 
multicomponent system replacement collisions also take 
place. Beside the contribution to the fast-neutron 
response associated with the production of electron-hole 
pairs, the TL output may also be affected by the produc- 
tion of lattice-ion displacements and lattice vibrations. A 
neutron- 

/erasing effect has been reported for LiF TLD phosphor 
(16 ). Fast neutrons seemed to have the effect of remo- 
ving electrons from their traps, thereby erasing some 
of the stored information, but Stanford and Johnson ( 3 ) 
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observed no significant indication of a neutron-erasing 

effect during their experiments. Goldstein, et al* (17) 

have found that for fast-neutron fluences as high as 
1 ^ 2 

10 y n/cm , the gamma-ray and fast neutron responses are 
additive, so that even though there are a sizable number 
of defects introduced, the TL output is not significantly 
affected. The fast-neutron- associated contribution to 
^LiF TL response was estimated by Simons and Yule (18) 
to be typically 10$ m a fast-reactor environment. Such 
contributions would obviously be different for different 
TLD materials. Puite (19) has reported that no permanent 
change m TL sensitivity for gamma-rays is observed after 
fast-neutron irradiations of upto 8 x 10 n/cm . 
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CHAPTER 3 


PRESENT APPROACH AND ANALYTICAL METHODS USED 

3.1 DESIGN OF EXPERIMENTS 

From the review of earlier gamma-heating 
measurements given m Section 2.2, it is seen that tempe- 
rature rises associated with radiation heating experiments 
m facilities such as low-power research reactors usually 
amount to very much less than a degree, necessitating 
the use of TLD's for deduction of energy deposition rates# 
Calorimetric methods involving temperature measurements 
are generally not applicable. The present work has 
involved the design of certain experiments to enable 
results of gamma-heating measurements made using TLD's 
to be directly assessed m the light of temperature mea- 
surements made using thermocouples. 

It was felt that for establishing the feasibility 
of such comparisons, a pure gamma - field provided by a 
strong cobalt - 60 source would be adequate. Accordingly, 
the FIPLY*f acility at B.A.R.C. (Section *f#1) which pro- 
vides a planar cobalt - 60 source of ~ 4*0 KCi strength, 
was chosen for the present study. At the irradiation 
position available ( 20 cms from the plane of the 

-Ip O mm i 

source), the gamma flux was # 10 cm sec , which is 
typical of that available with low-power reactors# 


♦Food Irradiation & Processing Laboratory 
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If one considers, for example, a 5 cm thick 
steel slab exposed to the above gamma flux, the order 
of magnitude of the steady - state temperature drop 
across the slab would be ^ 2 x 10**^ °C. This low value, 
of course, results from the high thermal conductivity 
of steel* The approach presently taken has been to use 
materials of low thermal conductivity m conjunction 
with materials of high thermal conductivity to obtain 
measurable temperature drops. Thus, to extend the above 
simple example of a 5 cm steel slab, consider this 
surrounded by 5 cm thick slabs of wood (a material often 
employed for neutron shielding). The steady - state 
temperature drop across the wood would be as high as 
^20 °C, if one assumes that all the heat generated m 
the steel is conducted through the wood. The gamma - 
heating m the high - conductivity steel thus gets 
reflected m an easily measurable temperature drop 
across the low-conductivity wood. 

In actual practice, du-=- to he at conduction 
through the steel m the Y and Z directions (taking X 
as the thickness direction), this temperature rise 
would not be ob earned . However, if insulating material 
is provided m the X and Z directions, so that most of 
the heat is m fact conducted through the wood, a 
measurable temperature drop would certainly be regis- 
tered. In the present work, 4 - - 8 cm of thermocole 
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was used for providing the necessary insulation for the 
experimental assemblies m the Y, Z directions (Sec. 

4.2) . 

A comparison of results for temperature diffe- 
rences deduced from dose measurements with actual thermo- 
couple measurements would involve an accurate heat - 
conduction calculation m the X, Y, Z direction. The 
3 - dimensional code, HEATING, was used for the purpose 
and this is described m Section 3*3* 

For obtaining energy deposition (heat genera- 
tion) rates at different locations m a given experimen- 
tal assembly using TLD’s, it is essential to estimate the 
gamma - energy spectrum at the various points (Sec. 2.3) » 
This was presently done by conducting one-dimensional 
photon transport calculations using the ANISN code 
(Sec. 3.2). 

While carrying out the temperature measurements, 
relatively long irradiations would be necessary for 
achieving steady - state conditions. With the FIPLY 
gamma - flux indicated earlier TLD irradiations, on the 
other hand, would have to be for typically ^ 2 mins, 
to avoid saturation in the response. For such short 
irradiations the effects of transient dose, received 
while bringing the source into position and while remo- 
ving it, must be accurately corrected for. Thus, in 
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the present case, the FIPLY source was located, under 
18 ft of water when not m use. For carrying out an 
irradiation the source was raised into position, an 
operation that took about 45 seconds during which signi- 
ficant doses were recorded by the TLD's. Similarly, at 
the end of the irradiation, a transient dose was recor- 
ded while lowering the source. Since gamma - heating 
rates had to be assessed m terms of watts per sec. of 
steady exposure, it was essential to carry out two sets 
of TLD measurements (for, say, 1 and 2 minutes of steady 
exposure) - so that the effects of the transients could 
be corrected for. 

3.2 ANISN CODE 

ANISN is a one-'Umensional, discrete-ordinates 
transport code and solves the one-dimensional, energy - 
dependent Boltzmann transport equation with anisotropic 
scattering for slab, cylindrical and spherical geometries 
(20, 21).ANISN solves forward or adjoint, homogeneous 
or inhomogeneous, neutronic or photonic problems. The 
inhomogeneous problems may have a fixed- volume, dis- 
tributed source or a specified angular - dependent 
shell source at ahy mesh interval. Fissions may be 
included for a subcntical system. Vacuum, reflec- 
tive, periodic, white or albedo boundary conditions 
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may be specified. Gross- sections may be input from a 
library tape and/or from cards. Fixed distributed 
sources or shell sources may be input from cards and/ 
or from tape. The code also includes space-point 
scaling to accelerate the flux solution on inner itera- 
tions . 

The discrete-ordinates, or Carlson's (22) 

S , method using a diamond - difference solution tech- 
nique is employed. The solution in the code will 
approach the exact solution of the Boltzmann equation 
with increasing orders of approximations as the space, 
angle and energy meshes approach differential size. 

The input data for the ANISN code are divided 
into the following five data sets; 

(A) Overall problem parameters 

(B) Cross-section data 

(C) Fixed- source data 

(D) Flux or fission guess data 

(E) Remainder of data. 

A discussion of (B), the cross-section data, 
is pertinent here. ANISN expects a table of cross- 
sections for each energy group, g, for each material in 
the following format : 



Position 


Cross-section type 


1 

Activity 



2 

* 

* 

u activity 

* 

■ 

1 Optional 

f 

I 

* 

IHT-2 

0 

'^absorption 

\ 

S 

j 

1 


XHT-1 

rr f ission 



IHT 

total 



IHT+1 

,r g , + NUS* •* 

g 


* 

« 

• 

• 


Up scatter 

0 

IHS-1 

°~g + 1 * g 



IHS 

•* g 



IHS+1 

0 

^g-1 * g 

V 


Downscatte 

0 

0 

IHM 

D 

* 

°g - KDS* 

i 

g _ 



Activities may Tee computed by zone and interval. 

This option provides a rapid and economical way to 
obtain data such as reaction rates, dose rates, cap- 
ture or other desired response rates. For the type 
of photonic calculations presently conducted (pair - 
production events negligible), 
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Absorption ~ £ (^he photoelectric cross-sections) 


<7 


total 


r + - r (the sura of photoelectric and 
Gompbon cross-sections) 


These values can be obtained from standard handbooks 
(23) at different gamma-energies. 


For the scattering cross-sections, with only 
down scatter possible, . t (g 1 > g) values are 

obtained from the Klein - Nishma equation (2*f). The 
equation is 

d KN r 2 ~2 

— — — = -5- 0 - f 1 + c< (1 - COS 0)1 . f 1 + cos 2 0 

dA ^ u J L 

2 2 

+ "- f +^'(l ' g -" S ~cos ~ e ' ) j cm2 /electron-sterdian 

(3.D 

E 

where oi~ y - FT T , r = 2.82 x 10" J cm (classical radius 
°' 511 0 of the electron) 

and 00s e = 1 + — f 1 1 - -• ■ -’ - W E being the initial 

B 0 E i o 

energy and E* the final energy after scattering through 
angle 0. 


Now, d-GL = 2 IT sm 0 d 0 . 

Putting (1 - cos 9) = z, so that sin 9 d 9 = dz, and 
integrating Equation (3*1) from z = 0 to = (1 - cos 9) 
we get 
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KN = r^- x 7,914. x 0.603 X 10 


+ •^*5’ ( 1 +C(Zj ) 

or 1 


-2 z£ 

A 

1 


2o( - c< 2 
1 


2o( ( 1 +.Xz,)^ (1 +o(zJ 


( + -4 ) - ln (1 +o(z 1 ).( ^ar + % - ~ )7 

or cx 3 1 a; 3 x °< - 1 

(3.2) 


cm 


Consider, for example, a ga.uma~energy group structure 

with 1 J+ - 1.1 Mev as the first group, 1.1 - 0.9 Mev 

as the second and so on (Section 5»2) . In this case 

<j: « is got from Equation (3-2) by taking E = 1.25 Mev 

1^1 o 

and E_ 1 =1.1 Mev. Eor <73 2 , E© = 1.25 Mev and 

E. =0.9 Mev. The other <7— values appropriate to 

the given group structure are obtained m a similar manner. 


3.3 HEATING CODE 

The HEATING program is a generalized heat- 
conduction code with reasonably simplified input require- 
ments (25) 26). The program establishes a network of 
nodal points throughout the specified configuration 
alongwith their relative locations and boundary conditions 
to enable the performing of nodal heat balances. The 
conductances from a given point to each of its neigh- 
bour points are calculated, and the program then compu- 
tes the temperature distribution for transient and/or 
steady-state heat-conduction problems m one, two or 
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three dimensions using either cartesian or cylindrical 
polar co-ordinates. 

The initial temperature pattern may he speci- 
fied as a function of position. Boundary temperatures 
may be time -dependent and heat generation rates may be 
time - and/or position - dependent. HEATING- offers nine 
possible geometries, which are members of cylindrical 
polar or rectangular co-ordinate systems, and three 
boundary types which are numbered 0, 1 or 2 corresponding 
to the following : 

0 - same effect as a zero-numbered boundary 

1 - convective heat transfer type 

2 - controlled surface temperature type. 

The spatial description of the problem is given m the 
form of mesh pomts/lmes m each direction upto a maxi- 
mum of 500 . 

The calculations arc based on the assumption 
that the heat transfer within the material follows 
Fourier's Law of heat conduction. In cartesian co- 
ordinates, the conduction equation for steady-state 
conditions is 

k t K (x > y ’ z) + % t K n z) §. 

+ dz [ K (Xs z) dzl* Q (x ’ y ’ ^ = 0 (3#3) 



where T is the temperature, K is the position- depen- 
dent thermal conductivity and Q is the position - 
dependent volumetric heat generation rate. 

The boundary condition is 

= f (x, y, z) (3 •*+- ) 

OR - K ( fE ) b = h [ T b - f (x, y, z)] (3.5) 

Consider a specified nodal point j which is thermally 

connected with K adjacent nodal points. By a heat 

balance on the node we get 

✓ 

M 

Z A (T - T ) + Q = 0 (3.6) 

J Jm J J 

Similar equations are written for all the nodes. The 

value of T is obtained by iteration. Assume that 
J 

n **“1 

the (n + 1)th iteration gives the correct value T^ ' 

and satisfies equation (3.6). The nth iteration value 

T 3 ? gives a residual value E when substituted m 
J 

equation (3.6), i.e. 

M 

1 It, (T_ - t“) + Q = B (3.7) 

1-1 J x x J J 

M ri4.-| 

Also, X K (T, - T^ +1 ) 
i=1 J 1 x J 

From Equations (3.7) and (3. 


+ q 3 = ° 


(3.8) 
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Rewriting Eqn, (3*7), we get 


M 


r 

1=1 


K 

D i 


(I x - 


r>n> 

'1 


Tt 1 ) + Q. 


^ rj.n+1 



M 

1=1 


oh 


( 3 . 10 ) 


To increase the convergence rate, one uses T = T^ for 
i > 3 and T l = T^ + ^ for l < 3 . By using Equation 
(3.10) the temperature imbalance is brought to minimum 
at all nodal points. An over- relaxation of the imbalance 
can be done by multiplying the left hand side of 
Equation (3*10) by a parameter /3 ,1.6. 


M 


M 


& [5 oh ‘h - + - *?> s, oh 

(3.11) 


Re-arranging the above equation M 

I K I + Q. 

T? +1 = (1 ~/3) + /3 ( - ) (3.12) 


O '" 1 _ i ' '3 

M 

JL A 

i=1 3 • L 


Use of Equation (3*12) is known as the extrapolated 
Liebmann method. /3 must be m the range 1 .0 - 2.0. The 
suggested optimum value of /? is 


0 = 2 



T n+2 _ T n+1 
T n+f _ T n 



1/2 


(3-13) 


y-j »pi .X.’l 

Equation (3.12) is iterated, replacing T by T , 

until the maximum absolute value of (T^ - T 1 ? 4 "^ )/ ( T 1 ? 4 ^ ) 

J J J 

is less than the specified convergence criterion, 

m lL 

say 10 . 



CHAPTER 4 


EXPERIMENTAL particulars 

The various experimental facilities used m the 
present work are described m this chapter. 

4.1 FI PLY FACILITY 

The cobalt-60 source employed for irradiation 
of the experimental assemblies was that of the FIPLY 
facility at the Bhabha Atomic Research Centre, Trombay . 
At the time of experimentation, the strength of the 
source was tsa 40 KCi. A schematic diagram of the FIPLY 
facility is shown in Fig. 4.1, The source consists of 
38 pencils m two racks as shown in Fig. 4.2 (a). Each 
pencil consists of a cluster of 5 cobalt-60 slugs (27), 
encapsulated m stainless steel as shown m Fig. 4.2 (b) . 
On either side of the source, irradiation channels of 
34 cm width and 26 cm height are provided (Fig. 4.1). 

The source room has staggered concrete shielding to 
minimize direct streaming and is separated from the 
control room by an interlocking door. Plate 1 gives 
a view of the FIPLY source room showing a "thermocouple 
assembly” (Sec. 4.2) being positioned m the irradia- 
tion channel. Plate 2 is a view of the control console 
m the FIPLY control room. 
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4.2 PLATES AND ASSEMBLIES 

Plates (300 mm length, 150 mm width and of 
various thicknesses) of steel, aluminium, lead and wood 
(teak) were fabricated to a specified surface finish of 
1.6 microns to ensure good thermal contact. Plate 3 
shows the various types of experimental plates. From 
these plates, different assemblies were made with ther- 
mocole insulation on each side, of 40 mm thickness m 
the Y - direction and 80 mm m the Z - direction. The 
various dimensions were largely dictated by the space 
available m the FIPLY irradiation channel* 

To enable dose measurements to be carried out, 
each plate m the "TLD assembly" was provided with nine 
16 mm - diameter holes (5 mm deep for steel and lead, 

8 mm for aluminium and 16 mm for wood). TLD - carrying 
blanks of the appropriate material and dimensions could 
be fitted in, these holes. Plate 4 gives a view of the 
various types of blanks used. Plate 5 shows the partly- 
assembled aluminium TLD assembly with thermocole 
insulation. 

For carrying out temperature measurements, it 
was felt that for ease of experimentation, a separate 
set of assemblies be used. These were constructed 
from thicker plates than employed for the TLD assemblies 
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so as to minimise any finite effects due to air gaps. 
Two of the plates of each material (excluding wood) 
were provided with alumel-chromel thermocouples, spot- 
welded on the surface at the centre with the leads 


running along 2 mm deep slots. Table 4.1 summarises 

j 

the configuration for each type of experimental assembly 


presently investigated. 

TABLE 4.1 The Various Assemblies 


No ♦ As s embly Overall Constituent Plates Constituent Plates 

Thickness for TLD Assemblies for Thermocouple 

Assemblies 

1 . Steel 160 mm 5 mm A1 *5 mm A1 

( 25 + 25 )mm W (25+25)mm W 

(6+6+6+6-+6+10+10)mm S ( $10+20+20 *)mm S 

(25+25)mm W ( 25+26 )mm W 

5mm AL 5* ram A1 


2 . Aluminium 


3. Lead 


210 mm 5 AL *5 mm A1 

(25+25)mm W (25+25)mm W 

( 10+1 0+1 0+1 0+20+20+20 )mm AL ( *20+20 +20+20+20 *)mm A1 
(25+25) mm W (25+ 2 ?)mm W 

5 mm AL 5* mm AL 


90 mm 5 mm AL 

(25+25) mm W 
(6+6+6+6+6) mm L 
(25+25) mm W 
5 mm A1 


+5 mm AL 
(25+25) mm W 

(*10+20*) mm L 
( 2 5+4 5.) mm W 
5* mm AL 


4. Combination-1 220 mm 5 mm AL 

(25+25)mm W 
(10+10)mm S 
(25+25)mm W 
( 10+1 0+1 0+1 0)mm AL 
(25+25)mm W 
5mm A1 


*5 mm A1 
( 25 + 25 )mm W 
(*20*)mm S 
( 25+25 )mm W 
(*20+20*)mm AL 
(25+25)mm W 
5*mm AL 


5. Combination^ 220 


mm 5mm AL 

(25+25)mm W 
( 10+1 0+1 0+1 0)mm A1 
( 25 + 25 )mm W 
(10+10)mm S 
(25+25)mm W 
5mm AL 


*5 mm AL 
(25+2 5 )mm W * 
(*20+20*) mm AL 
( 25+25 )mm W 


(*20*) mm S 
( 25+25 )mm W 
5* mm AL 
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In Table 4-.1 ^ * ' indicates a thermocouple 
mounted on the plate and Al, W, S and L denote 
Aluminium, Wood, Steel and Lead respectively. 

4-. 3 TLD'S AND TLD CALIBRATION FACILITY 

Pure CaFg and LiF powders prepared by the Health 
Physics Division, B.A.R.C., were the TLD's used. The 
parallel use of both these types of TLD was meant to 
provide a check on any possible systematic errors m 

-i 

calculating f - factors. 

If a dosimeter is to be used for absolute measu- 
rements, it must be calibrated m a known radiation field. 
During the calibration, as during the dose measurements, 
the TLD must be enclosed in a capsule of suitable thickness 
to establish electronic equilibrium (2). For cobalt-60 
radiations (mean energy, 1.25 Mev) , the amount of "build- 
up" material required is about 0.^5/^ cm., where -f is 
the density of the capsule material m gm/cm^. The build- 
up required for the various types of blanks used is com- 
pared m Table b.2 with the build-up actually provided. 

TABLE b.2 Build-up required and provided In TLD Blanks 



Build-un required 

Build-uu provided 

Steel 

0,57 mm 

1 .0 mm 

Aluminium 

1 .70 mm 

2.5 mm 

Lead 

0.4-0 mm 

1 .5 mm 

Wood 

5*60 mm 

6 .5 mm 
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The Elecbronics Division, BARC, has a standard 
point cobalt~60 source which was presently used for the 
calibration of the TLD‘s» This consists of a cylindrical 
(6 mm diameter, 6 mm height) capsule sealed m aluminium, 
the overall dimensions being 8 mm diameter and 9 mm height 
(28), The source is supported 4 ft. above ground level 
and located at the centre of a hall of dimensions 30 ft x 
30 ft x 30 ft to minimise scattering. The TLD calibration 
samples were placed 30 cm. away from the source at the 
same level and irradiated for a known time. The total 
exposure received by each sample was obtained in roentgens 
from the known source strength. The energy absorbed per 
unit mass by the TLD samples was then calculated m rads 
from 


E- 


’D 


0.88 + X c * 


^ 'm 


* f * e~' UX rads 


(4.1) 


•f air 

where 0,88 is the conversion factor for roentgens to rads. 
X c is gamma exposure m roentgens after correction 
for source decay. 


0 

( 




-”~) m is the mass energy absorption coefficient 
of capsule material at 1,2j?Mev. 

~~) is the mass energy absorption coefficient 
"T 9.1 X* 

of air at 1.25 Mev. 

: is the factor calculated by Equation (2.11) 
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*** 

and e is a correction for build-up material attenua- 

tion, x being the provided build-up thickness 
and ujl the linear attenuation coefficient « 

k.k TLD READER 

Several high-quality commercial TLD readers are 
available. Figure k.3 gives a schematic diagram for the 
system presently used. For carrying out a reading, the 
TLD powder was placed on the Kanthal strip which could be 
heated uniformly m a light-tight enclosure by a constant- 
current power supply. The TL emission produced was sensed 
by the photomultiplier tube, the signal being amplified 
by a D.C. amplifier and recorded on a Digilog recorder. 
Plate 6 gives a general view of the TLD reader set-up, 
while Plate 7 is a close-up view of the TLD - reader photo- 
multiplier tube with the TLD powder seen uniformly spread 
out on the Kanthal strip . 

The response of a TLD is obtained from the reader 
m tho form of a glow curve, i»e. TL mtensity-vs-tempe- 
rature/time. For talcing the reading, either the area 
under the glow curve or the maximum peak height may be 
taken as proportional to the dose. The latter parameter 
was presently used. Figure 4- A shows typical glow' curves 
obtained for CaF 2 and LiF in the course of the present 
experiments . 
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if. 5 TEMPERATURE MEASUREMENT SET-UP 

The temperature measurement set-up used m 
conjunction with each thermocouple assembly (Sec. 4.2) 
consisted of a potentiometer voltmeter bridge (PVB), a 
selector switch, and an ice bath (Plate 8). The PVB 
(e.s.i. Model 300) measured the potential difference 
developed across eoch thermocouple (4-0 microvolt/°C) to 
an accuracy of + 1 microvolt. The selector switch selec- 
ted the particular thermocouple being read. Ice used for 
the ice bath was prepared from distilled x^ater and the 
cold thermocouple ends dipped in mercury contained m 
separate test tubes inserted m the ice. Figure 4.5 shows 
the thermocouple connections for measurement of (a) absolute 
temperatures and (b) temperature differences. 
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FIG 43 SCHEMATIC DIAGRAM OF TLD READER 
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CHAPTER 5 


RESULTS AND THEIR DISCUSSION 
5.1 TLD'S USED 

Natural CaF^ and LiF TLD's m the form of 
powder were used for the present work. These were pre- 
pared Indigenously by the Health Physics Division, 

B.A.R.C » , the size of the grams for these powders being 
between 80 - 200 mesh. The effective densities of the 
CaE 2 and- powders were determined as 2.66 gm/cm^ and 
1.95 gm/cm , respectively. Standard annealing procedures 
were adopted before using the TLD's for dose measurements. 
Essentially the pre-ir radiation annealing involved heating 
of the TLD powder m an oven at 400 °C for one hour, 
followed by storage at 80 C for 24- hours. Fading tenden- 
cies of the TLD's were rendered negligible by the 
annealing . 

;.2 CALCULATION OF f g , f ± and f 

The deduction of energy absorbed by a material 

from dose measurements made using TLD's is dependent 

*1 

on the calculation of an appropriate ^ - factor for use 
in Equation (2.1) (Section 2.3). Equation (2.11) expre- 
sses f in terms of the small-cavity and large-cavity 

values f and f , respectively, for a given situation. 

S -L * 
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To calculate f g (T^ ), the electron- flux- 

dependent small -cavity values have to be averaged over 

the three primary-electron-spectrum components produced 

by the monoenergetic gamma-ray of energy (In the 

60 

present case of Co gammas, only photoelectric and 
Comp ton-interactions need be considered). Electron mass 
stopping powers to enable the deduction of f (Tv ) 

S 

from Equations (2.3) and (2.4) have been tabulated by 
Berger and Seltzer (29, 30). These references were 
presently used to provide the appropriate ( 18 ) mass 
stopping powers for steel, aluminium, lead and LiE. For 
CaF 2 , for which data was not directly available, stopping 
powers for aluminium and argon were suitably averaged 
since the effective atomic number (16*3) for CaF^ falls 
between these two elements. The stopping power for 
wood was obtained assuming an average weight composition 
of 50 % C, bb% 0 and 6 % H (31). 

For the evaluation of (T^ ) from Equation 

(2.6), tabulated values of mass energy absorption coeffi- 
cients for steel, aluminium and lead were taken from 
Reference (23). Values for wood were obtained by assu- 
ming the earlier-stated C/O/H composition. Reference 
(32) gave mass energy absorption coefficients for CaF 2 
and LiF. 
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Values of f g and f ± for CaF^steel, CaF^Al, 
GaF ^lead, CaF^wood, LiF/steel, LiF/Al, LiF/le a d and 
LiF/ wood were calculated for twelve different gajnma- 
energy values. These corresponded to the mean energies 
for the twelve photon-energy groups presently considered, 
m the range 0 .05-1 .40 Mev, for calculating the energy 
spectrum of the Co gamma-rays at different locations. 
Table 5*1 gives the gamma-energy group structure used, 
while Table 5*2 exemplifies the calculation of f g , f^ 
values for the case of CaF^steel. 

TABLE 5*1 Gamma- Energy Group Structure 


Energy Q Energy Range 11 Mean Energy 
Group jj (Mev) ( (Mev ) 


1 

1.4 - 1.1 

1.25 

2 

1.1 - 0.9 

1.00 

3 

0.9 - 0.7 

0.80 

4 

0 

* 

-vj 

1 

o 

O 

VJl 

0.60 

5 

0.5 - 0.3 

0.40 

6 

0.3 - 0.2 

0.25 

7 

0.2 - 0.15 

0.175 

8 

o.i5~ o.io 

0.125 

9 

0.10- 0.08 

0.090 

10 

0.08- 0.07 

0.075 

11 

0.07- 0.06 

0.065 

12 

0.06- 0.05 

0.055 
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The weighting factor d m Eqn. (2.11) was 
calculated by Eqn. (2.9), the effective mass absorption 
coefficient /3 being given by Equation ( 2 * 7 ) • The value 
of g used m Eqn. (2.9) was calculated for the LiF, 

Ca F 2 IbD's as exemplified by the following procedure 
for CaF 2 ; 


Average weight of CaF 2 powder 
obtained from dispenser, 

Effective density of CaF 2 powder, 
Volume of CaF 2 m blank, 

Diameter of hole in blank, 

Height of CaF 2 powder m blank, 


w ss 0.054- gm 
y = 2.66 gm/cm^ 

V = = 0.0203 crn^ 

i 

d = 0.3 mm (0.2 mm for 

lead blanks) 

h s -2- = 0 . 287 cm 


(0.64-6 cm for lead 
blanks) 


Average path length of electrons, g 


4- v 
s 


4- (% f d 2 h) 

( ‘| d 2 + "it a h) 


s 0.197 cm (0.173 for 
lead blanks) 

- 0.525 gm/cm 2 (0,460 

gm/cm 2 for lead blanks) 

The g value for LiF powder was similarly obtained as 
0.405 gm/cm 2 (0.281 gm/cm 2 for lead blanks). 

With the values of , g known, d was calcu- 
lated for each of the twelve gamma- energy values and 
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f (T v ) evaluated from Eqn. (2.11 ) for all the various 
combinations. The results are typified by Table 5.3 wmch 
gives *the values of y for CaF^steel , Curves giving 

the variation of f~^T^TJ with gamma-energy for the CaF^ 

TLD surrounded by steel, aluminium, lead and wood are 
shown m Fig, 5*1 , while Fig. 5.2 gives similar curves 
for the LiF TLD. 

TABLE 5*3 Calculation of y -/ 4 — r for CaFVsteel 
Combination 1 Uv ' * 


Energy $ 
Group 8 

/? 

1 1 r 

i a 5 

f 

s 

mm 


tfiU 

1/f (% ) 

1 

12.2 

0.156 

1.057 

0 *84^4* 

1 .068 

1.067 

0.937 

2 

16.8 

0.113 

1.058 

0.887 

1.061 

1 .061 

0.943 

3 

23.3 

0.082 

1.057 

0.918 

1.059 

1.059 

0.945 


35.7 

0.053 

1 .061 

0.947 

1 .031 

1.033 

0.968 

5 

66.1 

0.029 

1.067 

0.971 

0.977 

O. 98 O 

1 .021 

6 

138 

0.014 

1.073 

0 .986 

0.749 

0.75^ 

1.327 

7 

242 

0.008 

1.077 

0.992 

0.543 

0.547 

1 .829 

8 

471 

0.004 

1 .076 

0.996 

0.376 

0.379 

2.640 

9 

1000 

0.002 

1 .081 

0.998 

0 .314 

0.316 

3.167 

10 

1480 

0.001 

1.082 

0.999 

0.300 

0.301 

3.324 

11 

2290 

0.001 

1.079 

0.999 

0.298 

0.299 

3.348 

12 

4100 

0,000 

1.079 

1 .000 

0.299 

0.299 

3.340 


Notei The values of f are seen to be relatively insensitive 


to f s' 
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5.3 CALCULATION OF GAMMA- SPECTRA 

The various ejcperimental assemblies (Section 4,2) 
consisted of a number of individual plates of different 
materials. Gamma spectra at the centre of each plate 
were calculated using the ANISitf code (Section 3.2), S^, P Q 
approximations being presently employed. 


For modelling purposes, the FTPLY source (Fig. 

4.2) was taken as a distributed source of gamma-rays 

belonging to bhe first energy group, viz. 1.1 - 1.4 Mev. 

. 60 ,., 

(corresponding to the mean energy of 1.25 Mev for Go 
gamma’s). The homogenisation of the source region for 
obtaining appropriate photon cross-sections was carried 
out as follows s 


Area of 5 cobalt slugs of 0.635 cm. dia - 
Area of S.S. cladding — 
Total area = 
Area of 2.35 cm square 


1 . 58 cm 2 
2.02 cm 2 
3.60 cm 2 
5.52 cm 2 


Effective density of homogenised source q.60 A _ * 0 

region = x ~ • ? *' £ sm/w 

Photon cross-sections (absorption and total) 
for the source region (with the above effective density) , 
as well as for steel, aluminium, lead and wood (of 
measured densities) were obtained from Reference (23). 

Group scattering cross-sections for the different mate- 
rials were calculated using Eqn. (3.2) (Section 3.2). 
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Cross- sections for wood were calculated assuming the 
earlier mentioned weight composition of 5 0 % C, 0 and 

0> H, Table 5 .4 exemplifies the results for the group- 
to-group scattering cross-sections, the calculated values 
for steel being shown there. 

ANISN- calculated gamma- spectra at characteristic 
depths for the five different experimental assemblies 
(Sec. b. 2) are shown plotted m Figs. 5*3 - 5*7* 

A separate ANISN calculation was carried out to 

estimate the spectrum of gamma-rays emerging from the 

source region m the absence of any assembly. This is 

shown plotted m Fig. 5.7(a). The fraction of unscattered 

(Group 1) gamma-rays m this spectrum was 60 %, This 

is of the same order of magnitude as the results reported 

by Visser and Mouton (33), on the basis of experimental 

6o 

spectrum measurements of gamma-rays from a large Co 

source. Certain numerical checks were also carried out 

for the ANISN calculations, e.g., use of an Sg- approximation 

(instead of St) was found to have insignificant effect 

on the computed spectra. 

With the gamma- spectrum estimated at the centre 

of the various plates of each assembly, flux-weighted 

4 - values for each location were calculated using Eqn. 
f 

(2.12). Tables 5.5 - 5.9 give the jr - values calculated 
for various depths m the five different assemblies. 



TABLE 5.4- Gamma Scattering (Group-to-Group) Cross-Sections for Iron 
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Note? L.G, denotes the lower-energy group, for the scattering. 
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TABLE 5*5 Flux- weighted values of j in Steel 
Assembly for CaF 2 and LiF TLD'S 


Assembly 3 

Depth Ccm) 0 

Material | 

GaF 2 

if 

if 

LiF 

1.70 

Wood 

0.950 


1 .061 

4 .20 

Wood 

0.944 


1 .060 

5.35 

Steel 

1 .219 


1.584 

5.95 

Steel 

1 .214 


1 .524 

6.55 

Steel 

1 .212 


1.506 

7.15 

Steel 

1 .222 


1.520 

7.75 

Steel 

1 .226 


1*532 

8.45 

Steel 

1 .225 


1.527 

9.4-5 

Steel 

1 .236 


1.572 

11.7 

Wood 

0.930 


1.057 

14* • 2 

Wood 

0.926 


1 .051 



TABLE 5*6 Flux-weighted Values of ^ in Aluminium 
Assembly for CaF^ and LiF TLD'S 


Assembly 
Deoth (cm) 

| Material 

nr 

a 


nr 

5 

LiF 

1.70 

Wood 


0.892 


1.047 

4.20 

Wood 


0 .864 


1.042 

5.60 

A1 


0.841 


1.179 

6 .60 

A1 


0.828 


1 .201 

7.60 

A1 


0,816 


1 .219 

8.60 

A1 


0.806 


1.238 

10.60 

A1 


0.795 


1.260 

12.60 

A1 


0.7^3 


1.281 

14 .60 

A1 


0.777 


1.306 

16.70 

Wood 


0.724 


1.011 

19.20 

Wood 


0.732 


1.012 



TABLE 5*7 Flux-weighted Values of |r m Lead 
Assembly for CaF 2 and TLD'S 


Assembly 
Depth (cm) 

nr 

$„ 

Material | 

CaF 2 

nr 

a 

LiF 

1.70 


Wood 

0.987 


1,067 

4,20 


Wood 

0.992 


1.067 

5-35 


Lead 

1 .866 


1.683 

5.95 


Lead 

1.616 


1.432 

6.55 


Lead 

1.553 


1.368 

7.15 


Lead 

1.501 


1,326 

7.75 


Lead 

1.555 


1.375 

9.70 


WoGd 

1 .054 


1 .079 

12.20 


Wood 

1 .050 


1 .078 
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TABLE 5.8 Flux- weighted Values of ^ m Combination- 1 
for CaF 2 


Assembly 5 Material i CaJ 0 

Perth (cm) 8 8 £ 


Wood 

Wood 


1.70 

4.20 

5.^5 

6.45 

8.70 

11.20 

12.60 

13.60 

14.60 

15.60 

17.70 

20.20 


Steel 

Steel 

Wood 

Wood 

XL 

A1 

A1 

XL 

Wood 

Wood 


0.952 

0.948 

1.211 

1 .261 

0.880 

0.855 

0.836 

0.826 

0.820 

0.815 

0.774 

0.784 


1 \ f 

cei^o- j ' > ufijy 

54889 

Jke* No. * 





TABLE 5.9 Flux-weighted Values of jr m Combination- 2 
for CaE 2 


Assembly 
Depth (cm) 

| Material 

~r 

5 

CaF 2 

1.70 

Wood 


0.894 

4.20 

Wood 


0.868 

5.60 

Al 


0.844 

6 ,60 

A1 


0.833 

7.60 

Al 


0.822 

8 .60 

Al 


0.814 

10.70 

Wood 


0.794 

13*20 

Vfood 


0.813 

14.45 

Steel 


1.332 

15.45 

Steel 


1 .294 

17.70 

Wood 


0.884 

20,20 

Wood 


0.884 
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It is interesting to note that even in the 

case of a poorly matched combination, such as LiF/lead 

* 

(Fig. 5*2), the effective jr - values are not too far 
from unity. This is due to the filtering out of the 
low-energy gammas at greater depths within the high-Z 
material. This results m a harder spectrum for lead 
than for steel or aluminium (Figs. 5*3 - 5*5) and conse- 
quently m not too high a sensitivity on low-energy gammas 
for which j - variations are large. 

5.4 DOSE MEASUREMENT EXPERIMENTS 

Calibration of the TLD's was carried out before 

every assembly irradiation by filling four blanks each 

of steel» aluminium, lead and wood (as appropriate) with 

CaFgj and the same number of blanks with LiF . Half the 

blanks were placed at 30 cm, and the rest at 15 cm, from 

6o 

the standard Co calibration source discussed m Section 
4*3. Irradiation of these blanks was carried out for a 
known period of time ('-'4 days) to give doses m the same 
range as that acquired by TLD' s during the assembly irra- 
diations (10^ - 10^ rads). The doses absorbed by the 
calibration TLD’s were calculated by Eqn. (4.1), and 
their response obtained from the TLD reader under the same 
conditions as employed for the assembly TLD’ s. Typical 
dose/response curves obtained m the calibration process 



5 ^ 


are shown in Fig. 5,8, corrections for LiF sup ralme arity 
effects being thus explicitly obtained m uhe range of 
interest . 

It should be mentioned that for normalising the 
response of the TLD's irradiated m the different plates 
of each assembly, the calibration TLD's used were those 
irradiated m blanks of the corresponding material. This 
was to minimise the relative error m the j - corrections 
that had to be applied for deducing heat-generation rates 
m the plates. 

The FIPLY irradiation procedure for dose measure- 
ments in the various assemblies (Table 4-.1) was essentially 
the same as described below for the case of the steel 
assembly. 

Five wood blanks filled with TLD powder were 
placed m five different holes (the central and the four 
corner positions) m the first wood plate. The remaining 
holes m the plate carried wood blanks without TLD's. 
Similarly another five blanks with TLD's were located in 
the second wood plate, at the central and the four middle 
positions. The first steel plate carried five steel 
blanks containing TLD powder m the same positions as 
for the first wood plate, the second steel plate having 
five in the second set of positions. In this way* each 
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of the plates forming the assembly carried five TLD- 
contaimng blanks of the appropriate material, the set 
of positions alternating between the two mentioned above. 

The assembly was irradiated for 1 minute of 
steady exposure (with the source m the fully-up position)* 
After the irradiation, the TLD-contaimng blanks were taken 
out, and replaced by a second set of blanks containing fresh 
powder. A second irradiation of the assembly for 2 minutes 
of steady exposure was then carried out* 

The TLD response at each location corresponding 
to 1 minute of steady exposure had to be corrected for 
source-transient effects (Section 3*1 ) } and this was 
achieved by interpreting the various TLD readings as 
follows % 

If S denotes the corrected steady-exposure res- 
ponse (min~ 1 ) at a particular location, x the response 
due to source transients, X^ the reading obtained from 
the TLD used m the first irradiation and X £ that from 
the second, 

X^=S+x (5*1) 

X 2 = 2S + x (5*2) 

Equations (5.1) and (5.2) were used to eliminate x at 
each position, giving 


S 


(X 2 - X.,) 


(5.3) 
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From the appropriate calibration curve, the 
above response was expressed m rads/min for the TLD. 

The average result obtained for the five TLD locations 
m each plate gave the energy absorbed per minute for 
the TLD material m the plate. In the first few experi- 
ments, it was found that the average value was practically 
the same as that at the central location and, subsequently, 
only the central locations were employed. 

The heat generation m each plate was obtained 
m watts/cm^ of plate material by multiplying the above 
TLD dose by the appropriate -j: - factor (Section 5*2) and 
a conversion constant K given by 

~ Joule mm * _gm _ *rr5 * JL * ~ _ IQ — jp 

K “ ^dT * cm 3 - 10 60 f - 7 

(5*4) 

0 

being the density of the plate material m gm/cm . 

Table 5*16 gives the heat generation rates 
measured m the plates of the steel, aluminium and lead 
assemblies, using both CaF 2 and LiF . The results are shown 
graphically m Figs. 5*9 - 5*1 1» smooth curves having 
been drawn through the experimental points. These curves 
are shown normalised to unity at the centre of the assem- 
blies in Figs. 5.12 - 5.14, to enable comparison with 
gamma-heating results calculated by ANISN. The latter 
were obtained as a "reaction rate" from the ANISN - 



calculated gamma- spectra using "cross-sections" for the 
ith group given by 
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Z i = < yr 1 >i E i (5-0 

where ( — vr ) is the mass energy absorption coefficient 

T 1 

of the material for gammas of energy E , the mean energy 
of the ith group. 

TABLE 5*10 Heat-generation Rates (x 10“^) m Watts/cm^ 
Measured with CaFg and LiE m Steel, 
ALuminium and Lead Assemblies 


Steel Assembly jj 

Alluminium Ass 

embly | 

Lead Assembly 

Plate^j 

i CaP 2 8 

EH 

Plates | 

E8 

U 

11 

LOP | 

Plates 

Em 

LiF 

w i 

0.54 

0.41 

¥ 

o.4o 


0.42 

W 1 

0.28 

0.33 

4 

0,50 

0.38 

4 

0.38 


o.4o 

w 2 

0.26 

0.28 

4 

4. 2 

5.3 

A 1 

1.8 


1.6 

h 

9.5 

6.5 

s 2 

3.9 

3.7 

A 

1.5 


1.3 

L 2 

4.3 

4 • o 

S 3 

3.1 

3.0 

h 

1.1 


1.2 

L 3 

1.9 

2.2 


2.6 

2.8 


1.2 


1.1 

¥ 

1 #4 

1.8 

S 5 

2.1 

1.8 

a 5 

0.9 


0.9 

l 5 

1.3 

1.2 

H 

1.6 

1 .2 


0.6 


0.6 

W 3 

0.05 

0.05 

S 7 

1.2 

0.9 

Ay 

o.5 


0 .4 

\ 

0.05 

0.05 

w 3 

0.09 

0.05 

4 

0.12 


0.10 




\ 

0.08 

0.05 


0.11 


0.09 





(Notation for plates i W = wood, S = steel, A = aluminium, 

L = lead) . 
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It is seen from Table 5*10 that the heat genera- 
tion rates obtained from CaF 2 and LiF measurements are m 
reasonable agreement with each other, except for a few 
locations. While some of the differences may be attributed 
to random errors m the TLD measurements, significant con- 
tribution to the overall error may be expected from the 

•j 

fact that the j - factors have been deduced from infinite- 
slab geometry calculation of the gamma-spectra. This is 
illustrated, for example, by the heat generation values 
obtained from CaF 2 , LiF m the first lead plate of the 

lead assembly. With the gamma-spectrum being soft for the 

*1 

lead assembly at this location (Fig. 5*5), the j - factor 
for LiF is subject bo greater uncertainty than that for 
CaFg (Figs. 5.1, 5.2). Similar arguments, though to a 
smaller degree, would justify preference of the CaF 2 results 
over the LiF - measured values for steel and aluminium also. 
Accordingly, it was the CaFo ** measured heat-generation 
profiles of Figs. 5*9 - 5*11 that were used for providing 
the input to the HEATING code for the computation of tempe- 
rature drops. 

Heat generation rates m the various plates for 
Combinations 1 and 2 were measured only with CaF 2 for 
the above reasons. The results are tabulated m Table 
5.11 and shown graphically m Figs. 5*1 5 and 5*18. Also 
shown m the figures are the ANISN-c ale Hated heat- 
generation rates normalised to the experimental values at 
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the centre of the steel region for Combination- 1 , and. at 
the centre of the aluminium region for Combination-^ . 

It is seen from these figures, as also from Figs. 5,12- 
5.1 5 +, for the other assemblies, that the ANlSN-calculated 
heat-generation profiles are significantly different from 
the experimental profiles. This is only to be expected 
from the infinite-slab nature of the ABISN calculations. 


TABLE 5*11 Heat Generation Rates (x 10*"^) m 
Watts/cnr Measured with CaF 2 in 
Combination 1 and 2. 


Combination- 1 

;| Combination- 2 

Plates 

n 

5 Q 

| Plates | 

Q 

w 1 

0.31 

b 

0.4-3 

W 2 

0.30 

w 2 

0.32 

S 1 

3*9 

A i 

1.5 

S 2 

2.9 

A 2 

1 .4- 

W 3 

0.17 

a 3 

1 .2 


0.15 

\ 

1 .1 


0.30 

Vh 

0.20 

Ap 

0.29 

\ 

0.10 

b 

0,32 

S 1 

2.4 

\ 

0.19 

S 2 

1 .8 

w 5 

0.08 


0 .08 

W 6 

0.07 

W 6 

0.05 


(Notation for plates i W = wood ? S = steel, 

A = aluminium) 
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5*5 TEMPERATURE DROPS DEDUCED FROM DOSE MEASUREMENTS; 

Heat-generation rates in watts/ cm 3 obtained, 
from CaF^ measurements m each assembly were fed as 
part of the input data to the HEATING code. The density 
of each material was measured, and corresponding values 
of thermal conductivity were obtained from standard 
handbooks (3*+, 35} 36). Table 5*12 shows the density 
and thermal-conductivity values used for the various mate- 
rials . 

TABLE 5«t12 Density and Thermal Conductivity Values 
Used for the Various Materials 


Material jj 

Density 
gm/ cm3 

Q Thermal cond, 

5 watts/ cm °C 

Steel 

8.0 

0.52 (*) 

Aluminium 

2.8 

1.21 (*) 

Lead 

11.3 

0.35 (*) 

Wood 

0.72 

1 , 63 x 10“ 3 (**) 

Thermo cole 

0.017 

3.37X10" 4 " (***) 

* Ref. (3k) 

** Ref. 

(35) Ref. (36) 


Convective heat-transfer boundary conditions, 
with h = 11 x 10“ 3 watts/ cm 2 -°C and T aip = 30 °C, were 
used m the HEATING code. The pair of h, T alp values 
was obtained on the basis of consistency with measured 
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surface temperatures for the steel assembly.* Table 5*13 
summarizes the HEATING-calculated, steady-state T -values 
for the various assemblies, the figure below explaining 
the notation used for <£> T . 

(a) Steel, Aluminium, Lead, Assemblies 


FI PLY 
Source 


~y 






Wood 




-T -T 
12 1 


S t e el/ Al/ Lead 


A VV T 3 


Wood 






T. 


T, 


T. 


T, 


(b) Combinations 1 and 2 


PLY 

ource 


"Wood 


!<£* * 


Sueel/Al 


^T 0 =T ~T_ 
2 2 3 


Wood 




AL/Steel 


V/ood 


3 3 4- 


•1 


T, 


T. 


i 

Ti. 


<2\=T 4 -'i ? *Tj= t j- T 6 

I 

y 


3 


TABLE 5*13 Temperature Drops (°C) Obtained by HEATING Code 
Through Dose Measurements m Various Assemblies 


Assembly 

~r 

5 

h T 1 

Q 

5 

AI 2 

r 

_L 

£ T 0 

3 

<■ 

0<jO« 


Steel 


15.3 


0.02 


15.3 

**■ 


Aluminium 


11 .4- 


0.02 


11 .4- 


- 

Lead 


16.2 


0 03 


16.2 

- 

- 

Combination- 1 


14- .4 


0.003 


6,6 

o,oo4- 

8.0 

Combmation-2 


11 .8 


0.002 


1.1 

0.01 

10.3 


^Numerical studies vith H-UTIHG showed relatively little sensitivity 
of the calculated ji T - ’-lies th’ i-v^t i- 
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The HEATING results of TabD e 5*13 fully justify 
the design considerations for the present experiments 
(Sec. 3*1)j i.e., the temperature drops across the low- 
thermal-conductivity material, wood, are of a magnitude 
that can be easily compared with thermocouple measurements 
while those across steel, aluminium or lead are not, 

5.6 THERMOCOUPLE EXPERIMENTS s 

Each thermocouple assembly (Table 4.1) was irra- 
diated m FIPLY continuously for about 30 hours to achieve 
steady-state conditions. Extreme care was taken in assemb- 
ling the various plates prior to the irradiation to avoid 
any possibility of air gaps between them. Another precau- 
tion taken was to tape the front and rear surfaces of the 
assembly to avoid air circulation in any finite gaps between 
the thermocole insulation and the plates. Thermocouples 
readings were taken at 30 -mmute intervals during the 
irradiation following the procedure discussed m Section 
4*5* Both absolute temperatures and El T - values were 
obtained. Fig. 5.1? shows the ^ T-, transients for the 
steel, aluminium and lead assemblies, while Figs. 5»l8, 

5.19 give the 4^, A T^, T^ transients for Combina- 
tion 1 and 2, respectively. It is seen that steady-state 
conditions were achieved in all the experiments. 

The &. T-value s m Figs. 5»17**5.19 are seen to 
not pass through the origin. This was only to be expected 
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considering that the assembly was typically at <-'•30 °C, 
before being introduced into the air-cooled FI PLY irra- 
diation channel with an air-blast temperature of <-'25 °C. 
An experiment was conducted with a thermocouple assembly 
introduced into the FIPLY channel at time t = 0, the 
air-cooling being started without the source being brought 
out. An initial A T^ value of 5 °G was recorded for this 
experiment and this was found to drop to less than a 
degree within 8 hours as shown m Table 5«l4 ♦ Any initial 
AT effects would thus be quite insignificant after, 
say 20 hours. 

TABLE 5.14 Variation of Initial with Time 


Time m Hrs 

T 

X 

.AT.J in °G 

0 


5.0 

4 


1.9 

6 


1 .4 

8 


0.9 


Table 5»l5 summarises the steady-state A - 
values obtained during the irradiation of the various 
thermocouple assemblies. 
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TABLE 5*15 Temperature Drops (°C) Measured by 

Thermocouples in the Various Assemblies 


Assembly | 

* T 1 



mm 

Ed 


at 4 

T I — 

S 

Steel 

16.0 


0.20 


16,3 




Aluminium 

11 .9 


0.30 


11.9 


- 

- 

Lead 

17.0 


0.12 


17*2 


- 

- 

Combination- 1 

15.0 


0.03 


7*1 


0.10 

8.6 

Combmation-2 

12.5 


0.05 


2.7 


0.03 

10.5 


5.7 COMPARE SI ON OF DEDUCED AND MEASURED A 1 -VALUES 

From Tables 5.13 and 5*1 5? it is seen that the 
agreement between dose-deduced and measured temperature drops 
across the wood is quite satisfactory for all the assemblies, 
except for the AT^-results for Combination- 2. The latter 
discrepancy is explained by the following fact. VJhile m most 
of the other cases, HEATING-calculated results are more or 
less linearly dependent on the gamma-heating measurements for 
individual regions, the A T^- results for Combmation-2 depends 
on the difference of two nearly equal quantities, viz* the 
heat generated m the aluminium and steel regions for this 
assembly (Table 5*11)* 

Comparison of Tables 5*13 snd 5*1 5 does indi- 
cate a systematic trend, viz. that the theraiocouple-measured 
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temperature drops are about 5$ greater than those obtained 
from HEATING. The overall agreement is nevertheless quite 
satisfactory , considering the various sources of error 
discussed in the following section. 

5.8 DISCUSSION OE ERRORS 

Tables 5*16 and 5*17 summarize the sources of 
systematic and random error in temperature drops deduced 
from the gamma-heating measurements and those measured 
directly using thermocouples. The indicated estimates of 
the individual errors are based on order-of -magnitude 
considerations for a typical situation, e.g. error in 
HEATING-deduced A T-values have been assumed to be directly 
related to errors in measured heat— generation rates (not 
Valid for the A T^- results for Combmation-2) . 

TABLE 5.16 Sources of Error m Gamma-heating-deduced 
A T-V alues 


Source of Error 


Error 


(a) Systematic 


_ " -1 

1. Estimation of 4 - factors -f 2% 

2. Calibration of TLD's + 2% 

3. Thermal-conductivity data for wood and 

thermocole + 3/° 

HEATING code (boundary conditions etc.) + 2 % 

(b) Random (per experiment) 

1 . TLD Dispensing and Reading ± 3% 

2. Effect of Source Transients + 10 

3. Positioning of Assembly + 1% 



TABLE 5.17 Sources of Error in Thermocouple-measured 
AT-Values 


Source of Error 

Error 

(a) Systematic 


1. Effect of air gaps 

2. Reaching of steady-state 

3. Effect of fluctuations m outside 
temperature 

4 . The rmo c oupl e-End H e at mg 

5. Thermocouple Calibration 

+ 2% 

- ~\% 

± 1 % 

+ 

+ 1 \% 

(b) Random (per experiment) 


1 • Thermocouple Readout 

2. Positioning of Assembly 

negligibl 

± 


The overall (r.m.s) errors in the comparison of 
deduced and measured temperature drops (Section 5*7) have 
thus been presently estimated to be of the order of + 5$ 
systematic and + 3$ random. 
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CHAPTER 6 


CONCLUSIONS AND SCOPE FOR FURTHER WORK 
6.1 CONCLUSIONS 

The use of TLD’s for carrying oat gamma-heating 
measurements m several reactor materials has been assessed 
through temperature measurements across specially-designed 
assemblies. Thus, m effect, a direct comparison of 
dosimetric and calorimetric techniques for gauuna-heating 
evaluation has been shown possible. The materials consi- 
dered presently were steel, aluminium, lead and wood. 

The accuracy of gamma-heating determinations by 

the TLD technique is limited by uncertainties m the gamma- 

spectrum-dependent relationship between energy deposited m 

the dosimeter and that m the surrounding material (the 
"1 

- factor). The better the matching in effective atomic 
number between the dosimeter and surrounding material, the 
less is the uncertainty. Cobalt-60 gamma spectra within 
the experimental assemblies m the present study were esti- 
mated from infinite- slab-geometry calculations with the 

-i 

ANISN code. Adequacy of the calculated £ - factors was, 
however, established by obtaining reasonable agreement 
m energy-deposition rates deduced using LiF and CaFg 
TLD f s (with effective atomic numbers of 8.2 and 16.3, 
respectively) . 
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With the TLD-measured heat-generation rates fed 
as input , steady— state temperature drops across the low- 
conductivity material, wood, m the various assemblies 
were deduced from the three-dimensional heat-conduction 
code, HEATING. These were compared with actual thermocoupl 
measurements. The agreement obtained was within about 5 % 
m most cases, i.e. quite consistent with errors of + 5% 
systematic and + 3% random, as assessed for a typical 
experiment* 

6.2 SCOPE FOR FURTHER WORK 

The present study has been carried out using a 
large cobalt-60 source for irradiation. It would be 
instructive to carry out experiments of the present type 
m a reactor environment, so that 

(i) a more characteristic gamma- ray spectrum is 

considered for assessing the adequacy of gamma- 
heating measurements with ILD's and 
(ii) the accuracy of corrections for the neutron sensi 

i 

tivity of the TLD‘s may also be assessed. 

Accurate estimation of gamma spectra (parti- 
cularly the low-energy component) would be essential 
for evaluating appropriate -j: - factors. Theoretical 
est im ates could be supported by experimental measurements, 
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e.g. the comparison of the TLD response, in both the 
reactor and. the calibration facility, with ion-current 
chamber measurements (3). 

for minimal neutron sensitivity, a higher effec- 
tive atomic number for the TLD would probably be better 
because of the lower average energy of the recoil ions 

produced In elastic scattering. The preference of Cal 2 

7 

over Lif for gamma-heating measurements thus becomes even 
more justified m a reactor environment. (^LiF would, of 
course, nob be usable because of its high neutron absorption 
cross-section) , 
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